Insights into the Epidemiology, Phylodynamics, and Evolutionary Changes of Lineage GI-7 Infectious Bronchitis Virus

Infectious bronchitis virus (IBV) is distributed worldwide and causes significant losses in the poultry industry. In recent decades, lineages GI-19 and GI-7 have become the most prevalent IBV strains in China. However, the molecular evolution and phylodynamics of the lineage GI-7 IBV strains remain largely unknown. In this study, we identified 19 IBV strains from clinical samples from January 2021 to June 2022 in China, including 12 strains of GI-19, 3 strains of GI-7, and 1 strain each of GI-1, GI-9, GI-13, and GI-28. These results indicated that lineages GI-19 and GI-7 IBVs are still the most prevalent IBVs in China. Here, we investigated the evolution and transmission dynamics of lineage GI-7 IBVs. Our results revealed that the Taiwan province might be the origin of lineage GI-7 IBVs and that South China plays an important role in the spread of IBV. Furthermore, we found low codon usage bias of the S1 gene in lineage GI-7 IBVs. This allowed IBV to replicate in the host during evolution as a result of reduced competition, mainly driven by natural selection and mutational pressure, where the role of natural selection is more prominent. Collectively, our results reveal the genetic diversity and evolutionary dynamics of lineage GI-7 IBVs, which could assist in the prevention and control of viral infection.

Identification of Evolutionary Trajectories Shared across Human Betacoronaviruses.

Comparing the evolution of distantly related viruses can provide insights into common adaptive processes related to shared ecological niches. Phylogenetic approaches, coupled with other molecular evolution tools, can help identify mutations informative on adaptation, although the structural contextualization of these to functional sites of proteins may help gain insight into their biological properties. Two zoonotic betacoronaviruses capable of sustained human-to-human transmission have caused pandemics in recent times (SARS-CoV-1 and SARS-CoV-2), although a third virus (MERS-CoV) is responsible for sporadic outbreaks linked to animal infections. Moreover, two other betacoronaviruses have circulated endemically in humans for decades (HKU1 and OC43). To search for evidence of adaptive convergence between established and emerging betacoronaviruses capable of sustained human-to-human transmission (HKU1, OC43, SARS-CoV-1, and SARS-CoV-2), we developed a methodological pipeline to classify shared nonsynonymous mutations as putatively denoting homoplasy (repeated mutations that do not share direct common ancestry) or stepwise evolution (sequential mutations leading towards a novel genotype). In parallel, we look for evidence of positive selection and draw upon protein structure data to identify potential biological implications. We find 30 candidate mutations, from which 4 (codon sites 18121 [nsp14/residue 28], 21623 [spike/21], 21635 [spike/25], and 23948 [spike/796]; SARS-CoV-2 genome numbering) further display evolution under positive selection and proximity to functional protein regions. Our findings shed light on potential mechanisms underlying betacoronavirus adaptation to the human host and pinpoint common mutational pathways that may occur during establishment of human endemicity.

The Mutational Landscape of SARS-CoV-2.

Mutation research is crucial for detecting and treating SARS-CoV-2 and developing vaccines. Using over 5,300,000 sequences from SARS-CoV-2 genomes and custom Python programs, we analyzed the mutational landscape of SARS-CoV-2. Although almost every nucleotide in the SARS-CoV-2 genome has mutated at some time, the substantial differences in the frequency and regularity of mutations warrant further examination. C>U mutations are the most common. They are found in the largest number of variants, pangolin lineages, and countries, which indicates that they are a driving force behind the evolution of SARS-CoV-2. Not all SARS-CoV-2 genes have mutated in the same way. Fewer non-synonymous single nucleotide variations are found in genes that encode proteins with a critical role in virus replication than in genes with ancillary roles. Some genes, such as spike (S) and nucleocapsid (N), show more non-synonymous mutations than others. Although the prevalence of mutations in the target regions of COVID-19 diagnostic RT-qPCR tests is generally low, in some cases, such as for some primers that bind to the N gene, it is significant. Therefore, ongoing monitoring of SARS-CoV-2 mutations is crucial. The SARS-CoV-2 Mutation Portal provides access to a database of SARS-CoV-2 mutations.

Molecular detection of monkeypox and related viruses: challenges and opportunities.

The recent widespread emergence of monkeypox (mpox), a rare and endemic zoonotic disease by monkeypox virus (MPXV), has made global headlines. While transmissibility (R0 ≈ 0.58) and fatality rate (0-3%) are low, as it causes prolonged morbidity, the World Health Organization has declared monkeypox as a public health emergency of international concern. Thus, effective containment and disease management require quick and efficient detection of MPXV. In this bioinformatic overview, we summarize the numerous molecular tests available for MPXV, and discuss the diversity of genes and primers used in the polymerase chain reaction-based detection. Over 90 primer/probe sets are used for the detection of poxviruses. While hemagglutinin and A-type inclusion protein are the most common target genes, tumor necrosis factor receptor and complement binding protein genes are frequently used for distinguishing Clade I and Clade II of MPXV. Problems and possibilities in the detection of MPXV have been discussed.

Epidemiological and Molecular Evolution of the Hiv-1 Crf94: Birth of Crf132

Background: In 2018 we reported the emergence of the new HIV-1 recombinant CRF94-02BF2 involved in a large transmission cluster of 49 French MSM mostly infected in 2016-2017. This CRF94 raised concerns of enhanced virulence. Prevention actions were undertaken in the area and population affected. This study reported the molecular and epidemiological evolution of this CRF94 until June 2022. Method(s): In 2021-2022, French sequence databases were screened for patients infected with HIV-1 subtype CRF94 or similar strain. HIV subtyping was confirmed by phylogenetic analysis of genes encoding both protease and reverse transcriptase (1070bps), and integrase (696bps) using IQ-Tree. Five whole genomes, related but distinct from CRF94, were obtained with the DeepChek assay Whole Genome kits. Recombination breakpoints were estimated using RDP4 and SimPlot. Mann-Whitney and LogRank tests were used for statistical analyses to compare patients' characteristics. Result(s): In June 2022, 49 new HIV-1 sequences were collected: 14 clustered with the 49 previous CRF94, 32 formed a new cluster next to but distinct from CRF94, and 3 strains could not be classified. Analysis of 5 whole genomes from the new cluster revealed a new recombinant, the CRF132-94B, mainly consisting of CRF94 which recombined with subtype B in the POL and accessory genes. Vif gene changed from the F2 to the B subtype. Both CRF94 and 132 clusters involved >95% of MSM, mostly infected < 1 year before diagnosis. However, there were differences: 97% were diagnosed in 2013-2019 for CRF94 vs 90% in 2020-2022 for CRF132. At time of diagnosis, 33% of patients infected with CRF94 knew the Prep vs 95% for CRF132. In the cluster CRF94, patients were older (34 vs 30 years, p=0.02), had higher viral loads (5.42 vs 4.42 log10 copies/Ml;p< 0.001), a lower CD4 cell counts (358 vs 508 /mm3, p=0.002). On treatment, the patients with the CRF94 reached viremia < 50 copies/Ml significantly later than those infected with CRF132 (p=0.0002). The prevention activities targeting the CRF94 cluster could explained the few patients infected with this strain after 2018. The CRF132 is mainly located in another Paris region area, but no specific transmission place has been identified. Conclusion(s): After 2019, the CRF94 spread seems greatly slowed down but the very close CRF132-94B has given birth to a new highly active cluster in 2020- 2022, despite the COVID social-distancing and a strong knowledge of the Prep. CRF132 appears to be less virulent perhaps due to the Vif gene change. Identified breakpoints positions of the new HIV-1 CRF132-94B. GenBank accession numbers of the five references : ON901787 to ON901791.

Molecular and Structural Evolution of Porcine Epidemic Diarrhea Virus.

To analyze the evolutionary characteristics of the highly contagious porcine epidemic diarrhea virus (PEDV) at the molecular and structural levels, we analyzed the complete genomes of 647 strains retrieved from the GenBank database. The results showed that the spike (S) gene exhibited larger dS (synonymous substitutions per synonymous site) values than other PEDV genes. In the selective pressure analysis, eight amino acid (aa) sites of the S protein showed strong signals of positive selection, and seven of them were located on the surface of the S protein (S1 domain), suggesting a high selection pressure of S protein. Topologically, the S gene is more representative of the evolutionary relationship at the genome-wide level than are other genes. Structurally, the evolutionary pattern is highly S1 domain-related. The haplotype networks of the S gene showed that the strains are obviously clustered geographically in the lineages corresponding to genotypes GI and GII. The alignment analysis on representative strains of the main haplotypes revealed three distinguishable nucleic acid sites among those strains, suggesting a putative evolutionary mechanism in PEDV. These findings provide several new fundamental insights into the evolution of PEDV and guidance for developing effective prevention countermeasures against PEDV.

Adaptive Evolution of the Spike Protein in Coronaviruses.

Coronaviruses are single-stranded, positive-sense RNA viruses that can infect many mammal and avian species. The Spike (S) protein of coronaviruses binds to a receptor on the host cell surface to promote viral entry. The interactions between the S proteins of coronaviruses and receptors of host cells are extraordinarily complex, with coronaviruses from different genera being able to recognize the same receptor and coronaviruses from the same genus able to bind distinct receptors. As the coronavirus disease 2019 pandemic has developed, many changes in the S protein have been under positive selection by altering the receptor-binding affinity, reducing antibody neutralization activities, or affecting T-cell responses. It is intriguing to determine whether the selection pressure on the S gene differs between severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and other coronaviruses due to the host shift from nonhuman animals to humans. Here, we show that the S gene, particularly the S1 region, has experienced positive selection in both SARS-CoV-2 and other coronaviruses. Although the S1 N-terminal domain exhibits signals of positive selection in the pairwise comparisons in all four coronavirus genera, positive selection is primarily detected in the S1 C-terminal domain (the receptor-binding domain) in the ongoing evolution of SARS-CoV-2, possibly owing to the change in host settings and the widespread natural infection and SARS-CoV-2 vaccination in humans.

Molecular evolution and structural analyses of the spike glycoprotein from Brazilian SARS-CoV-2 genomes: the impact of selected mutations.

The COVID-19 pandemic caused by SARS-CoV-2 has reached by February 2022 more than 380 million cases and 5.5 million deaths worldwide since its beginning in late 2019, leading to enhanced concern in the scientific community and the general population. One of the most important pieces of this host-pathogen interaction is the spike protein, which binds to the hACE2 cell receptor, mediates the membrane fusion and is the major target of neutralizing antibodies against SARS-CoV-2. The multiple amino acid substitutions observed in this region, specially in RBD have enhanced the hACE2 binding affinity and led to several modifications in the mechanisms of SARS-CoV-2 pathogenesis, improving the viral fitness and/or promoting immune evasion, with potential impact in the vaccine development. In this work, we identified 48 sites under selective pressures, 17 of them with the strongest evidence by the HyPhy tests, including VOC related mutation sites 138, 142, 222, 262, 484, 681, and 845, among others. The coevolutionary analysis identified 28 sites found not to be conditionally independent, such as E484K-N501Y. The molecular dynamics and free energy estimates showed the structural stabilizing effect and the higher impact of E484K for enhanced binding affinity between the spike RBD and hACE2 in P.1 and P.2 lineages (specially with L452V). Structural changes were also identified in the hACE molecule when interacting with B.1.1.7 RDB. Despite some destabilizing substitutions, a stabilizing effect was identified for the majority of the positively selected mutations.Communicated by Ramaswamy H. Sarma.

Functional binding dynamics relevant to the evolution of zoonotic spillovers in endemic and emergent Betacoronavirus strains.

Comparative functional analysis of the dynamic interactions between various Betacoronavirus mutant strains and broadly utilized target proteins such as ACE2 and CD26, is crucial for a more complete understanding of zoonotic spillovers of viruses that cause diseases such as COVID-19. Here, we employ machine learning to replicated sets of nanosecond scale GPU accelerated molecular dynamics simulations to statistically compare and classify atom motions of these target proteins in both the presence and absence of different endemic and emergent strains of the viral receptor binding domain (RBD) of the S spike glycoprotein. A multi-agent classifier successfully identified functional binding dynamics that are evolutionarily conserved from bat CoV-HKU4 to human endemic/emergent strains. Conserved dynamics regions of ACE2 involve both the N-terminal helices, as well as a region of more transient dynamics encompassing residues K353, Q325 and a novel motif AAQPFLL 386-92 that appears to coordinate their dynamic interactions with the viral RBD at N501. We also demonstrate that the functional evolution of Betacoronavirus zoonotic spillovers involving ACE2 interaction dynamics are likely pre-adapted from two precise and stable binding sites involving the viral bat progenitor strain's interaction with CD26 at SAMLI 291-5 and SS 333-334. Our analyses further indicate that the human endemic strains hCoV-HKU1 and hCoV-OC43 have evolved more stable N-terminal helix interactions through enhancement of an interfacing loop region on the viral RBD, whereas the highly transmissible SARS-CoV-2 variants (B.1.1.7, B.1.351 and P.1) have evolved more stable viral binding via more focused interactions between the viral N501 and ACE2 K353 alone.Communicated by Ramaswamy H. Sarma.

Genomic characterization and molecular evolution of SARS-CoV-2 in Rio Grande do Sul State, Brazil.

SARS-CoV-2 is the virus responsible for the COVID-19 and has afflicted the world since the end of 2019. Different lineages have been discovered and the Gamma lineage, which started the second wave of infections, was first described in Brazil, one of the most affected countries by pandemic. Therefore, this study analyzed SARS-CoV-2 sequenced genomes from Esteio city in Rio Grande do Sul, Southern Brazil. We also comparatively analyzed genomes of the two first years of the pandemic from Rio Grande do Sul state for understanding their genomic and evolutionary patterns. The phylogenomic analysis showed monophyletic groups for Alpha, Gamma, Delta and Omicron, as well as for other circulating lineages in the state. Molecular evolutionary analysis identified several sites under adaptive selection in membrane and nucleocapsid proteins which could be related to a prevalent stabilizing effect on membrane protein structure, as well as majoritarily destabilizing effects on C-terminal nucleocapsid domain.

Boosters for adeno-associated virus vector (AAV) (r)evolution.

Adeno-associated virus (AAV) is one of the most exciting and most versatile templates for engineering of gene-delivery vectors for use in human gene therapy, owing to the existence of numerous naturally occurring capsid variants and their amenability to directed molecular evolution. As a result, the field has witnessed an explosion of novel "designer" AAV capsids and ensuing vectors over the last two decades, which have been isolated from comprehensive capsid libraries generated through technologies such as DNA shuffling or peptide display, and stratified under stringent positive and/or negative selection pressures. Here, we briefly highlight a panel of recent, innovative and transformative methodologies that we consider to have exceptional potential to advance directed AAV capsid evolution and to thereby accelerate AAV vector revolution. These avenues comprise original technologies for (i) barcoding and high-throughput screening of individual AAV variants or entire capsid libraries, (ii) selection of transduction-competent AAV vectors on the DNA level, (iii) enrichment of expression-competent AAV variants on the RNA level, as well as (iv) high-resolution stratification of focused AAV capsid libraries on the single-cell level. Together with other emerging AAV engineering stratagems, such as rational design or machine learning, these pioneering techniques promise to provide an urgently needed booster for AAV (r)evolution.

Molecular evolution of SARS-CoV-2 from December 2019 to August 2022.

BACKGROUND: Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) pandemic spread rapidly and this scenario is concerning worldwide, presenting more than 590 million coronavirus disease 2019 cases and 6,4 million deaths. The emergence of novel lineages carrying several mutations in the spike protein has raised additional public health concerns worldwide during the pandemic. AIM: The present study review and summarizes the temporal spreading and molecular evolution of SARS-CoV-2 clades and variants worldwide. The evaluation of these data is important for understanding the evolutionary histories of SARS-CoV-2 lineages, allowing us to identify the origins of each lineage of this virus responsible for one of the biggest pandemics in history. METHOD: A total of 2,897 SARS-CoV-2 whole-genome sequences (WGSs) with available information from the country and sampling date (Dec/2019 - Aug/2022), were obtained and were evaluated by Bayesian approach. RESULTS: The results demonstrated that the SARS-CoV-2 the time to the most recent common ancestor (tMRCA) in Asia was 2019-12-26 [Highest Posterior Density 95% (HPD95%): 2019-12-18; 2019-12-29), in Oceania 2020-01-24 (HPD95%: 2020-01-15; 2020-01-30), in Africa 2020-02-27 (HPD95%: 2020-02-21; 2020-03-04), in Europe 2020-02-27 (HPD95%: 2020-02-20; 2020-03-06), in North America 2020-03-12 (HPD95%: 2020-03-05; 2020-03-18) and in South America 2020-03-15 (HPD95%: 2020-03-09; 2020-03-28). Between Dec 2019 and June 2020, 11 clades were detected [20I (Alpha) and 19A, 19B, 20B, 20C, 20A, 20D, 20E (EU1), 20F, 20H (Beta)]. From July to Dec 2020 four clades were identified [20J (Gamma, V3), 21C (Epsilon), 21D (Eta), and 21G (Lambda)]. Between Jan and June 2021, three clades of the delta variant were detected (21A, 21I, and 21J). Between July and Dec 2021, two variants were detected, delta (21A, 21I, and 21J) and omicron (21K, 21L, 22B, and 22C). Between Jan and June 2022, the delta (21I and 21J) and omicron (21K, 21L, and 22A) variants were detected. Finally, between July and Aug 2022 three clades of omicron were detected (22B, 22C, and 22D). Clade 19A was first detected in the SARS-CoV-2 pandemic (Wuhan strain) with origin in 2019-12-16 (HPD95%: 2019-12-15; 2019-12-25); 20I (Alpha) in 2020-11-24 (HPD95%: 2020-11-15; 2021-12-02); 20H (Beta) in 2020-11-25 (HPD95%: 2020-11-13; 2020-11-29); 20J (Gamma) was 2020-12-21 (HPD95%: 2020-11-05; 2021-01-15); 21A (Delta) in 2020-09-20 (HPD95%: 2020-05-17; 2021-02-03); 21J (Delta) in 2021-02-26 (2020-11-02; 2021-04-24); 21M (Omicron) in 2021-01-25 (HPD95%: 2020-09-16; 2021-08-08); 21K (Omicron) in 2021-07-30 (HPD95%: 2021-05-30; 2021-10-19); 21L (Omicron) in 2021-10-03 (HPD95%: 2021-04-16; 2021-12-23); 22B (Omicron) in 2022-01-25 (HPD95%: 2022-01-10; 2022-02-05); 21L in 2021-12-20 (HPD95%: 2021-05-16; 2021-12-31). Currently, the omicron variant predominates worldwide, with the 21L clade branching into three (22A, 22B, and 22C). Phylogeographic data showed that Alpha variant originated in the UK, Beta in South Africa, Gamma in Brazil, Delta in India, Omicron in South Africa, Mu in Colombia, Epsilon in the USA, and Lambda in Peru. CONCLUSION: The COVID-19 pandemic has had a significant impact on global health worldwide and the present study provides an overview of the molecular evolution of SARS-CoV-2 lineage clades (from the Wuhan strain to the currently circulating lineages of the omicron). This article is protected by copyright. All rights reserved.

Molecular Evolution of SARS-CoV-2 during the COVID-19 Pandemic.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) produced diverse molecular variants during its recent expansion in humans that caused different transmissibility and severity of the associated disease as well as resistance to monoclonal antibodies and polyclonal sera, among other treatments. In order to understand the causes and consequences of the observed SARS-CoV-2 molecular diversity, a variety of recent studies investigated the molecular evolution of this virus during its expansion in humans. In general, this virus evolves with a moderate rate of evolution, in the order of 10-3-10-4 substitutions per site and per year, which presents continuous fluctuations over time. Despite its origin being frequently associated with recombination events between related coronaviruses, little evidence of recombination was detected, and it was mostly located in the spike coding region. Molecular adaptation is heterogeneous among SARS-CoV-2 genes. Although most of the genes evolved under purifying selection, several genes showed genetic signatures of diversifying selection, including a number of positively selected sites that affect proteins relevant for the virus replication. Here, we review current knowledge about the molecular evolution of SARS-CoV-2 in humans, including the emergence and establishment of variants of concern. We also clarify relationships between the nomenclatures of SARS-CoV-2 lineages. We conclude that the molecular evolution of this virus should be monitored over time for predicting relevant phenotypic consequences and designing future efficient treatments.

The evolution of the spike protein and hACE2 interface of SARS-CoV-2 omicron variants determined by hydrogen bond formation.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was first detected in December 2019. As of mid-2021, the delta variant was the primary type; however, in January 2022, the omicron (BA.1) variant rapidly spread and became the dominant type in the United States. In June 2022, its subvariants surpassed previous variants in different temporal and spatial situations. To investigate the high transmissibility of omicron variants, we assessed the complex of spike protein 1 receptor-binding domain (S1RBD) and human angiotensin-converting enzyme 2 (hACE2) from the Protein Data Bank (6m0j, 7a91, 7mjn, 7v80, 7v84, 7v8b, 7wbl and 7xo9) and directly mutated specific amino acids to simulate several variants, including variants of concern (alpha, beta, gamma, delta), variants of interest (delta plus, epsilon, lambda, mu, mu without R346K) and omicron variants (BA.1, BA.2, BA.2.12.1, BA.4, BA.5). Molecular dynamics (MD) simulations for 100 ns under physiological conditions were then performed. We found that the omicron S1RBD-hACE2 complexes become more compact with increases in hydrogen-bond interactions at the interface, which is related to the transmissibility of SARS-CoV-2. Moreover, the relaxation time of hydrogen bonds is relatively short among the omicron variants, which implies that the interface conformation alterations are fast. From the molecular perspective, PHE486 and TYR501 in omicron S1RBDs need to involve hydrogen bonds and hydrophobic interactions on the interface. Our study provides structural features of the dominant variants that explain the evolution trend and their increased contagiousness and could thus also shed light on future variant changes.

A Probabilistic Approach to Evaluate the Likelihood of Artificial Genetic Modification and Its Application to SARS-CoV-2 Omicron Variant

A method to find a probability that a given bias of mutations occur naturally is proposed to test whether a newly detected virus is a product of natural evolution or a product of non-natural process such as genetic manipulation. The probability is calculated based on the neutral theory of molecular evolution and binominal distribution of non-synonymous (N) and synonymous (S) mutations. Though most of the conventional analyses, including dN/dS analysis, assume that any kinds of point mutations from a nucleotide to another nucleotide occurs with the same probability, the proposed model takes into account the bias in mutations, where the equilibrium of mutations is considered to estimate the probability of each mutation. The proposed method is applied to evaluate whether the Omicron variant strain of SARS-CoV-2, whose spike protein includes 29 N mutations and only one S mutation, can emerge through natural evolution. The result of binomial test based on the proposed model shows that the bias of N/S mutations in the Omicron spike can occur with a probability of 2.0 × 10−3 or less. Even with the conventional model where the probabilities of any kinds of mutations are all equal, the strong N/S mutation bias in the Omicron spike can occur with a probability of 3.7 × 10−3, which means that the Omicron variant is highly likely a product of non-natural process including artifact. © 2022 Information Processing Society of Japan.

RNA structure-altering mutations underlying positive selection on Spike protein reveal novel putative signatures to trace crossing host-species barriers in Betacoronavirus.

Similar to other RNA viruses, the emergence of Betacoronavirus relies on cross-species viral transmission, which requires careful health surveillance monitoring of protein-coding information as well as genome-wide analysis. Although the evolutionary jump from natural reservoirs to humans may be mainly traced-back by studying the effect that hotspot mutations have on viral proteins, it is largely unexplored if other impacts might emerge on the structured RNA genome of Betacoronavirus. In this survey, the protein-coding and viral genome architecture were simultaneously studied to uncover novel insights into cross-species horizontal transmission events. We analysed 1,252,952 viral genomes of SARS-CoV, MERS-CoV, and SARS-CoV-2 distributed across the world in bats, intermediate animals, and humans to build a new landscape of changes in the RNA viral genome. Phylogenetic analyses suggest that bat viruses are the most closely related to the time of most recent common ancestor of Betacoronavirus, and missense mutations in viral proteins, mainly in the S protein S1 subunit: SARS-CoV (G > T; A577S); MERS-CoV (C > T; S746R and C > T; N762A); and SARS-CoV-2 (A > G; D614G) appear to have driven viral diversification. We also found that codon sites under positive selection on S protein overlap with non-compensatory mutations that disrupt secondary RNA structures in the RNA genome complement. These findings provide pivotal factors that might be underlying the eventual jumping the species barrier from bats to intermediate hosts. Lastly, we discovered that nearly half of the Betacoronavirus genomes carry highly conserved RNA structures, and more than 90% of these RNA structures show negative selection signals, suggesting essential functions in the biology of Betacoronavirus that have not been investigated to date. Further research is needed on negatively selected RNA structures to scan for emerging functions like the potential of coding virus-derived small RNAs and to develop new candidate antiviral therapeutic strategies.

A molecular evolution algorithm for ligand design in DOCK.

As a complement to virtual screening, de novo design of small molecules is an alternative approach for identifying potential drug candidates. Here, we present a new 3D genetic algorithm to evolve molecules through breeding, mutation, fitness pressure, and selection. The method, termed DOCK_GA, builds upon and leverages powerful sampling, scoring, and searching routines previously implemented into DOCK6. Three primary experiments were used during development: Single-molecule evolution evaluated three selection methods (elitism, tournament, and roulette), in four clinically relevant systems, in terms of mutation type and crossover success, chemical properties, ensemble diversity, and fitness convergence, among others. Large scale benchmarking assessed performance across 651 different protein-ligand systems. Ensemble-based evolution demonstrated using multiple inhibitors simultaneously to seed growth in a SARS-CoV-2 target. Key takeaways include: (1) The algorithm is robust as demonstrated by the successful evolution of molecules across a large diverse dataset. (2) Users have flexibility with regards to parent input, selection method, fitness function, and molecular descriptors. (3) The program is straightforward to run and only requires a single executable and input file at run-time. (4) The elitism selection method yields more tightly clustered molecules in terms of 2D/3D similarity, with more favorable fitness, followed by tournament and roulette.

Exceptional diversity and selection pressure on coronavirus host receptors in bats compared to other mammals.

Pandemics originating from non-human animals highlight the need to understand how natural hosts have evolved in response to emerging human pathogens and which groups may be susceptible to infection and/or potential reservoirs to mitigate public health and conservation concerns. Multiple zoonotic coronaviruses, such as severe acute respiratory syndrome-associated coronavirus (SARS-CoV), SARS-CoV-2 and Middle Eastern respiratory syndrome-associated coronavirus (MERS-CoV), are hypothesized to have evolved in bats. We investigate angiotensin-converting enzyme 2 (ACE2), the host protein bound by SARS-CoV and SARS-CoV-2, and dipeptidyl-peptidase 4 (DPP4 or CD26), the host protein bound by MERS-CoV, in the largest bat datasets to date. Both the ACE2 and DPP4 genes are under strong selection pressure in bats, more so than in other mammals, and in residues that contact viruses. Additionally, mammalian groups vary in their similarity to humans in residues that contact SARS-CoV, SARS-CoV-2 and MERS-CoV, and increased similarity to humans in binding residues is broadly predictive of susceptibility to SARS-CoV-2. This work augments our understanding of the relationship between coronaviruses and mammals, particularly bats, provides taxonomically diverse data for studies of how host proteins are bound by coronaviruses and can inform surveillance, conservation and public health efforts.

Evolutionary dynamics of the severe acute respiratory syndrome coronavirus 2 genomes.

The coronavirus disease 2019 (COVID-19) pandemic has caused immense losses in human lives and the global economy and posed significant challenges for global public health. As severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, has evolved, thousands of single nucleotide variants (SNVs) have been identified across the viral genome. The roles of individual SNVs in the zoonotic origin, evolution, and transmission of SARS-CoV-2 have become the focus of many studies. This review summarizes recent comparative genomic analyses of SARS-CoV-2 and related coronaviruses (SC2r-CoVs) found in non-human animals, including delineation of SARS-CoV-2 lineages based on characteristic SNVs. We also discuss the current understanding of receptor-binding domain (RBD) evolution and characteristic mutations in variants of concern (VOCs) of SARS-CoV-2, as well as possible co-evolution between RBD and its receptor, angiotensin-converting enzyme 2 (ACE2). We propose that the interplay between SARS-CoV-2 and host RNA editing mechanisms might have partially resulted in the bias in nucleotide changes during SARS-CoV-2 evolution. Finally, we outline some current challenges, including difficulty in deciphering the complicated relationship between viral pathogenicity and infectivity of different variants, and monitoring transmission of SARS-CoV-2 between humans and animals as the pandemic progresses.

Chirohepevirus from Bats: Insights into Hepatitis E Virus Diversity and Evolution

Homologs of the human hepatitis E virus (HEV) have been identified in more than a dozen animal species. Some of them have been evidenced to cross species barriers and infect humans. Zoonotic HEV infections cause chronic liver diseases as well as a broad range of extrahepatic manifestations, which increasingly become significant clinical problems. Bats comprise approximately one-fifth of all named mammal species and are unique in their distinct immune response to viral infection. Most importantly, they are natural reservoirs of several highly pathogenic viruses, which have induced severe human diseases. Since the first discovery of HEV-related viruses in bats in 2012, multiple genetically divergent HEV variants have been reported in a total of 12 bat species over the last decade, which markedly expanded the host range of the HEV family and shed light on the evolutionary origin of human HEV. Meanwhile, bat-borne HEV also raised critical public health concerns about its zoonotic potential. Bat HEV strains resemble genomic features but exhibit considerable heterogeneity. Due to the close evolutionary relationships, bat HEV altogether has been recently assigned to an independent genus, Chirohepevirus. This review focuses on the current state of bat HEV and provides novel insights into HEV genetic diversity and molecular evolution.